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Observation of Field-Induced Transverse Neel Ordering 
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Neutron elastic scattering experiments have been performed on the spin gap system TlCuCls 
in magnetic fields parallel to the 6-axis. The magnetic Bragg peaks which indicate the field- 
induced Neel ordering were observed for magnetic field higher than the gap field H s « 5.5 
T at Q = (h, 0, 1) with odd I in the a* — c* plane. The spin structure in the ordered phase 
was determined. The temperature and field dependence of the Bragg peak intensities and the 
phase boundary obtained were discussed in connection with a recent theory which describes the 
field-induced Neel ordering as a Bose-Einstein condensation of magnons. 
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The singlet ground state with the excitation gap (spin 
gap) is a notable realization of the macroscopic quan- 
tum effect in quantum spin systems. When a mag- 
netic field is applied in the spin gap system, the gap 
A is suppressed and closes completely at the gap field 
H s = A/<7/ib- For H > H g the system can undergo 
magnetic ordering due to three-dimensional (3D) inter- 
actions. Such field-induced magnetic ordering was stud- 
ied first for Cu(NC>3)2 ■ IH2O. 1 ) However, the magnetic 
properties near H = H g have not been investigated be- 
cause of the very low ordering temperature, the maxi- 
mum of which is 0.18 K. Recently, the study of field- 
induced magnetic ordering has been revived, because the 
spin gap has been found in many quantum spin systems. 
Field-induced 3D ordering has been observed in several 
quasi-one-dimensional spin gap systems. 2 ~ 6 ' 

This paper is concerned with field-induced magnetic 
ordering in TICUCI3. This compound has a monoclinic 
structure (space group P2i/c). 7 ) TlCuCl 3 contains pla- 
nar dimers of CU2CI6, in which Cu 2+ ions have spin-i 
These dimers are stacked on top of one another to form 
infinite double chains parallel to the crystallographic a- 
axis. These double chains are located at the corners and 
center of the unit cell in the b — c plane, and are sep- 
arated by Tl + ions. The magnetic ground state is the 
spin singlet with the excitation gap A/fee ~ 7.5 K. 8,9 ) 
The magnetic excitations in TICUCI3 were investigated 
by Oosawa et a/., 10 ' who found that the lowest excita- 
tion occurs at Q — (0, 0, 1) and its equivalent reciprocal 
points, as observed in KCUCI3. 11, 12 -* The origin of the 
gap is the strong antiferromagnetic interaction J = 5.26 



meV on the planar dimer CU2CL3 in the double chain. 
The neighboring dimers couple magnetically along the 
chain and in the (1,0,-2) plane. 

Our previous magnetic measurements revealed that 
TICUCI3 undergoes 3D magnetic ordering in magnetic 
fields higher than the gap field H g w 5.5 T. 9 ) The mag- 
netization exhibits a cusplike minimum at the ordering 
temperature Tr. The phase boundary on the tempera- 
ture vs field diagram is independent of the field direction 
when normalized by the g-i actor, and can be represented 
by the power law 



{3 /2) [H n (T)-H s ]kT<I> 



(1) 



with 4> = 2.2, where H^(T) is the transition field at 
temperature T. These features cannot be explained by 
the mean-field theory. 13,14 ) 

Nikuni et al. 15 ^ argued that field- induced magnetic or- 
dering in TICUCI3 can be represented as a Bose-Einstein 
condensation (BEC) of excited triplets (magnons). In 
the magnon BEC theory, the external field H and the 
total magnetization M are related to the chemical po- 
tential /1 and the number of magnons N, respectively, as 
(j. = gnB(H — H g ) and M = g/JsN. The observed tem- 
perature dependence of the magnetization and the field 
dependence of the ordering temperature 9 ^ were qualita- 
tively well described by the magnon BEC theory based 
on the Hartree-Fock (HF) approximation. 15 ) 

If the magnons undergo BEC at ordering vector Qo for 
H > H g , then the transverse spin components have long- 
range order, which is characterized by the same wave vec- 
tor Qo- When the magnetization is small, i.e., magnons 
are dilute, the transverse magnetization per site is 
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Tabic I. Atomic coordinates of TICUCI3. 





X 


y 


z 


Tl 


0.7780 


0.1700 


0.5529 


Cu 


0.2338 


0.0486 


0.1554 


Cl(l) 


0.2656 


0.1941 


0.2597 


Cl(2) 


0.6745 


-0.0063 


0.3177 


Cl(3) 


-0.1817 


0.0966 


-0.0353 



expressed by 



m± = 5MB \/n c /2 , 



where n c is the condensate density. Thus the condensate 
density n c is proportional to the magnetic Bragg peak 
intensity. 

In order to confirm the transverse spin ordering cor- 
responding to the off-diagonal long-range order of the 
BEC, and to investigate the temperature and field depen- 
dence of the transverse magnetization m± corresponding 
to the condensate density n c , we performed neutron elas- 
tic scattering experiments on TICuCb;. 

The details of sample preparation were reported in 
reference 9. Since the precise structural analysis for 
TICUCI3 has not been reported, we performed the single- 
crystal X-ray diffraction. The lattice parameters at 
room temperature were determined as a = 3.9815A, 
b = 14.1440A, c = 8.8904A and (i = 96.32°. The atomic 
coordinates obtained are shown in Table I. 

Neutron scattering experiments were carried out at 
the El spectrometer of the BER II Research Reactor 
of the Hahn-Meitner Institute with the vertical field cry- 
omagnet VM1. The incident neutron energy was fixed at 
Ei = 13.9 meV, and the horizontal collimation sequence 
was chosen as 40'-80'-40'. A single crystal of <~ 0.4 cm 3 
was used. The sample was mounted in the cryostat with 
its (0, 1, 0) cleavage plane parallel to the scattering plane, 
so that the reflections in the a* — c* plane were investi- 
gated. The sample was cooled to 0.05 K using a dilution 
refrigerator. An external magnetic field of up to 12 T was 
applied along the 6-axis. Lattice parameters a = 3.899 
A, c = 8.750 A and (3 = 95.55° were used at helium 
temperatures. 

Magnetic Bragg reflections were observed for H > 
H g « 5.5 T, which could be indexed by Q — (h, 0, 1) 
with odd I. Figure 1 shows the 8 — 26 scans for (1, 0, —3) 
reflection measured at H — and 12 T at 1.9 K. These 
reciprocal lattice points are equivalent to those for the 
lowest magnetic excitation at zero field. 10 ' This result 
indicates that the magnetic unit cell is the same as the 
chemical one. Ferromagnetic Bragg reflections due to 
the induced moment could not be detected in the present 
measurements. The integrated intensities of nine mag- 
netic reflections were measured at 1.9 K under an exter- 
nal field of 12 T. The results are summarized in Table II 
together with the calculated intensities. The intensities 
of nuclear Bragg reflections were also measured to deter- 
mine the magnitude of the magnetic moment per site. 
The nuclear peaks were observed at Q = (h, 0, 1) with 
even as expected from the space group P2i/c. How- 
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Fig. 1. 9 — 20 scans for (1, 0, —3) reflection measured at H = 
and 12 T at 1.9 K. 



Table II. Observed and calculated magnetic Bragg peak intensi- 
ties at T = 1.9 K and H = 12 T for H\\b. The intensities are 
normalized to the (0, 0, 1)m reflection. R is the reliability factor 



given by R = Y, h ,k,l l J cal 
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(h,k,l) 


-^obs 


-fcal 


(0, 0, 1) M 


1 


1 


(0, 0, 3) M 





0.014 


(0, 0, 5) M 


0.063 


0.107 


(1, 0, 1) M 


0.031 


0.009 


(1, 0, -1) M 


0.074 


0.127 


(1, 0, 3) M 


0.056 


0.024 


(1, 0, -3) M 


0.398 


0.380 


(2, 0, 1) M 


0.012 


0.013 


(2, 0, -1) M 


0.128 


0.102 


R 


0.12 





ever, very weak nuclear peaks were also observed for odd 
I. To refine the magnetic structure, we used the atomic 
coordinates shown in Table I and the nuclear scattering 
lengths &ti = 0.878, b Cu = 0.772 and bci = 0.958 with 
the unit of 10~ 12 cm. 16 ' The magnetic form factors of 
Cu 2+ were taken from reference 17. The extinction ef- 
fect was evaluated by comparing observed and calculated 
intensities for various nuclear Bragg reflections. 

Figure 2 shows the spin structure determined at 1.9 K 
and 12 T. The directions of spins are in the a — c plane 
which is perpendicular to the applied field. Spins on the 
same dimers represented by thick lines in Fig. 2 are an- 
tiparallel. Spins are arranged in parallel along a leg in the 
double chain, and make an angle of 39° with the a-axis. 
The spins on the same legs in the double chains located 
at the corner and the center of the unit cell in the b — c 
plane are antiparallel. Although scans in the reciprocal 
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Fig. 2. Spin structure in the ordered phase of TICUCI3. The ex- 
ternal field is applied along the fe-axis. The double chain located 
at the corner and the center of the chemical unit cell in the b — c 
plane are represented by solid and dashed lines, respectively. 



plane including the 6*-axis were not performed, the spin 
structure shown in Fig. 2 is uniquely determined within 
the present measurements. Thus the transverse magnetic 
ordering predicted by the theory 13-15 ) was confirmed. 

Comparing magnetic peak intensities with those of 
nuclear reflections, we evaluate the magnitude of the 
transverse magnetization at 1.9 K and 12 T as mj_ = 
gHs(S±) = 0.26(2) fiQ. Since the Neel temperature 
for H = 12 T is T N = 7.2 K, the magnitude of to_l 
at 1.9 K must be very close to that at zero tempera- 
ture (see Fig. 2). For T — > 0, the condensate density 
n c approximates the density of the magnon n, which is 
given by n = 2m/<?/XB with the magnetization m per 
site. The value of n at 12 T and 1.8 K is obtained as 
n w 0.028 with m « 0.03 and g = 2.06. 8 ' 9 ) Sub- 
stituting n c s» n 1=3 0.028 into eq. (2), m± « 0.24 ^b- 
This value agrees with m± = 0.26(2) obtained by the 
present measurement. 

Figure 3 shows the temperature dependence of the 
magnetic peak intensity for (1,0,-3) measured in var- 
ious magnetic fields. On the right abscissa, we show 
which corresponds to the condensate density n c = 
2 (mi/ g^B.) 2 ■ The Neel ordering can be detected at the 
temperatures indicated by the thick arrows. With in- 
creasing external field, the ordering temperature T^{H) 
and the intensity of the magnetic Bragg peak increase. 
The anomaly around T^(H) becomes smeared as the 
magnetic field approaches H g , which leads to the dif- 
ficulty in determining the Neel temperature from the 
temperature dependence of the Bragg peak intensity for 
H < 6.6 T. Since the sharp anomaly is observed at 
the transition field Hn(T) in the field dependence of 
the Bragg intensity (see Fig. 4), and the sharp phase 
transition was detected through the specific heat mea- 



Fig. 3. Temperature dependence of the Bragg peak intensity for 
(1, 0, —3) reflection in various magnetic fields. The square of 
the transverse magnetization per site is shown on the right 
abscissa. 

surement, 18 ) the smearing of the peak intensity around 
Tn(H) is attributed not to the intrinsic smearing of the 
ordering temperature due to additional effects such as 
the staggered inclination of the principal axis of the g- 
tensor and the Dzyaloshinsky-Moriya interaction, but to 
the diffuse scattering. The amount of diffuse scattering 
around T^(H) seems to be independent of peak intensity 
at T = 0. 

Figure 4 shows the field dependence of the Bragg inten- 
sity for (1, 0, —3) reflection measured at various tempera- 
tures. On the right abscissa, we show m\ . The intensity 
remains almost zero up to the transition field (T) , and 
then increases rapidly. It is clear that field-induced Neel 
ordering takes place at H^(T). At 0.2 K, the intensity 
is almost proportional to H — H g just above i? g , 
and is slightly concave at higher fields. With increasing 
temperature, the transition field Hjq(T) increases, and 
the slope just above H^(T) increases. 

The magnon BEC theory based on the HF approxi- 
mation 15 ) predicts that the magnon density n and the 
condensate density n c are both concave functions of 
pi = g^i B (H - H s ) for T^0. 15) The deviation from the 
linear function of /i is due to the noncondensate magnons, 
and is proportional to /i 3 / 2 . The field dependence of 
the Bragg intensity for sufficiently low temperatures is 
explainable by the magnon BEC theory. However, the 
magnetization for H > H g is the slightly convex func- 
tion of H, 8 ^ as observed in low-dimensional spin systems. 
This behavior disagrees with the magnon BEC theory. 

We summarize the phase transition temperature 
Tn(H) and field Hn(T) in Fig. 5, in which the data ob- 
tained from the previous magnetic measurements are also 
plotted. Because the transition points determined from 
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Fig. 4. Magnetic field dependence of the Bragg peak intensity 
for (1, 0, —3) reflection at various temperatures. The square of 
the transverse magnetization per site is shown on the right 
abscissa. 



the present neutron scattering experiments and the mag- 
netic measurements lie almost on the same line, they are 
consistent with each other. This supports the predic- 
tion of the magnon BEC theory that the magnetization 
exhibits a cusplikc minimum at the transition tempera- 
ture. 15) 

The phase boundary can be expressed by the power law 
of eq. (1). In the previous magnetic measurements, the 
lowest temperature was 1.8 K. The low-temperature data 
contribute greatly to the determination of the exponent 
({>. Therefore, adding the present data down to 0.05 K, we 
reevaluate the value of <j>. We fit eq. (1) to the data for 
T < 4 K, which is lower than half the gap temperature 
A/k B = 7.5 K. The best fit is obtained with H g = 5.54 
T and <j> — 2.0(1), which is slightly smaller than the 
previous value of <f> = 2.2. The exponent obtained, <p — 
2.0(1), is somewhat larger than the <j> = 1.5 predicted by 
the magnon BEC theory. The difference between these 
values may be attributed to the fluctuation effect which 
is disregarded in the HF approximation. 

In conclusion, we have presented the results of neutron 
elastic scattering experiments performed on the spin gap 
system TICUCI3 in high magnetic fields along the 6-axis. 
The field-induced transverse Neel ordering was clearly 
observed for H > H g w 5.5 T. The phase boundary 
can be described by the power law as predicted by the 
magnon BEC theory. However, the exponent obtained, 
(j) = 2.0, is somewhat larger than the predicted value of 
= 1.5. 

The authors would like to thank T. Nikuni and M. Os- 
hikawa for stimulating discussions. H. T., A. O. and T. 
K. express their gratitude to the Hahn-Meitner Institute 
for their kind hospitality during the experiment. 



Fig. 5. The phase boundary in TICUCI3 for H\\b determined from 
the results of the temperature variation (closed rectangles) and 
field variation (closed circles ) of the Bragg peak intensity. Open 
circles and rectangles denote the transition points determined 
from the previous magnetization measurements. The dashed line 
is a guide for the eyes. The solid line denotes the fit by eq. (1) 
with (g/2)H c (0) = 5.7(1) T and <f> = 2.0(1). 
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